Many-electron dynamics triggered by massively parallel ionization 
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Massively parallel ionization of many atoms in a cluster or bio-molecule is identified as new phe- 
nomenon of light-matter interaction which becomes feasible through short and intense FEL pulses. 
Almost simultaneously emitted from the illuminated target the photo-electrons can have such a 
high density that they interact substantially even after photoionization. This interaction results in 
a characteristic electron spectrum which can be interpreted as convolution of a mean-field electron 
dynamics and binary electron-electron collisions. We demonstrate that this universal spectrum can 
be obtained analytically by summing synthetic two-body Coulomb collision events. Moreover, we 
propose an experiment with hydrogen clusters to observe massively parallel ionization. 
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Modern light sources such as free-electrons lasers (FELs) 
[1, 2] couple large numbers of photons into clusters [3] 
or bio-molecules [4], or more generally, extended sys- 
tems. Within femtoseconds many electrons are released 
through single-photon absorption and the ions left be- 
hind form a deep background potential. In cases, where 
most electrons are trapped in this potential, one ob- 
serves a sharp transition from continuous equilibration of 
the photo-activated electrons [5-8] to a non-equilibrium 
plasma executing characteristic oscillations [9] as the 
pulse length falls below the relaxation time. If the elec- 
trons are activated with sufficient energy to escape the 
potential in large numbers, a similar transition occurs 
when the pulse length falls below a critical escape time 
enabling direct interaction and energy exchange among 
the electrons even after photo-ionization, as we will show 
in the following. The regime introduced and discussed 
here is the exact opposite to the previously investigated 
case of "non-interacting electrons" in sequential emission 
(also referred to as multi-step ionization [10]), which oc- 
curs for long pulses [7, 10, 11]. 

The high energy of the excited electronic system per- 
mits a treatment in terms of classical Coulomb dynamics 
of ions and electrons [3, 4, 12, 13]. This is a tremendous 
simplification and allows us to calculate the time evo- 
lution of this many-body system and the photo-electron 
spectrum (PES) which results from illumination of the 
cluster with an intense high-frequency laser pulse using 
classical molecular dynamics with photo-ionization rates 
for the atomic ionization within the atomic cluster [12]. 

We will interpret the PES in two different ways: In 
terms of global types of dynamics we will show that mas- 
sively parallel ionization can be thought of as a convo- 
lution of a mean-field component and a component typ- 
ical for binary collisions. In terms of detailed paths of 
electrons we will demonstrate that the electron spectrum 
can be reproduced extremely well by approximating each 
photo-electron's final energy through a sum of contribu- 
tions from synthetic binary collisions with each of the 
other electrons. 

This two-fold interpretation is facilitated by the con- 



cept of Coulomb complexes (CC), which we have intro- 
duced recently [7] . They capture the essentials of electron 
dynamics activated through multiple photo-ionization in 
a cluster or bio-molecule. To understand massively paral- 
lel ionization, we only need CCs in their simplest version: 
A single isolated photoionization event from the mother 
atom in the cluster leads to an excess energy of E* . UN 
photo-electrons are produced by a laser pulse of length T 
the ions left behind are assumed to form a smooth back- 
ground potential, see Fig. 1. It is Coulombic outside the 
(spherical) cluster of radius i?, V{r) — ~N/r for r > R, 
and parabolic within the cluster, 

V{r) = Vo [rV(3i?2) - l] , r < R (1) 

with the depth Vq = "^N/R. Activated electrons 
are propagated under this potential and their mutual 
Coulomb repulsion. The overall electronic dynamics 
of the CC is completely determined by the four pa- 
rameters {N,R,E*,T). Moreover, CCs are scale in- 
variant, i.e., the one-parameter manifold of CC(?7) = 
{(iV, r]R, ?r'^E*, r]-^/'^T) \ ?7>0} leads to the same scaled 
dynamics. On the one hand, this renders phenomena 
which can be described by CCs quite general and, on the 
other hand, facilitates to identify a parameter combina- 
tion which can be realized in an experiment. 




FIG. 1: Sketch of A'^-fold photo-ionization of atoms in a cluster 
of radius R, leading to an ionic back ground potential (blue) 
with total charge and depth Vo = 3N/2R, cf. Eq. (1). 
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In Fig. 2 we compare the PES from the full molecu- 
lar dynamics calculation (where electrons and ions move 
classically according to all Coulomb forces) with the one 
obtained using the CC with its static and smooth ionic 
background. Obviously, both results agree with each 
other quite well demonstrating that CCs are a realistic 
approximation for the present scenario. Motivated by the 
scaling property of the Coulomb complex, we rescale the 
energy by the depth of the ionic background potential 
Vo = ^N/R = ?7 of all activated electrons potential. In 
Fig. 2b the PES is plotted in terms of the scaled energy 
e ^'q-'^E = E/Vo. 

The form of the PES exhibits some resemblance to the 
mean-field result (dashed lines in Fig. 2) which can be 
obtained analytically from the CC: Assuming as before 
spherical geometry, the potential for an electron photo 
activated with excess energy E* at radial distance r' from 
the center is Vynf(r') = V{r') + Vcc{r'), where V is given 
in Eq. (1) and Vcc{r') — Nr'/r' is the repulsive potential 
of the charge Nr' — Nr'^/R^ created by all electrons 
within the sphere of radius r' . To escape from the CC 
the electron has to overcome T^^f and its final energy is 
therefore E = E* + Vaif{r') or in scaled units 
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With the radial electron distribution dP/dr' — 3r'^/R^ 
and Eq. (2) we get for dP/ds = [dP/dr'] [de/dr']-'^ 
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within the interval e* — \ < e < £* which is of length unity 
or Vb in unsealed energies. The width Vq of the PES gives 
an account of the depth of the potential and consequently 
of the charge to extension ratio of the Coulomb complex. 
In particular the full width at half maximum 



l^E=^-Va=^-N/R or 



(4) 



of the mean-field spectrum is quite similar to its coun- 
terpart in the full spectrum. In the regime of massively 
parallel ionization this result is very useful to estimate the 
number of photons absorbed if the cluster size is known, 
or vice versa, determine the size of the cluster illuminated 
if one can measure how many electrons (their number 
equals the number of photons N) have been released. 

Apart from the overall agreement one observes in Fig. 2 
that the accurate PES is substantially blurred at the 
boundaries compared to the mean-field spectrum. The 
broadening can be interpreted as the result of a convo- 
lution with a spectrum governed by binary collisions in- 
duced by a s/iort-range, singular potential, i.e., the exact 
opposite of mean-field dynamics, which is generated by 
smooth long-vsaige interaction. 
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For the sake of being specific we model the short-range 
interaction of two electrons with a distance r by a Yukawa 
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FIG. 2: Color online: Photo-electron spectrum for sudden 
massively parallel ionization, a) Ari47 with icosahedral ge- 
ometry. Two thirds of the atoms loose an electron from the 
3p level with an excess energy of i5* = 0.4keV. b) Coulomb 
complex with = 100 and an excess energy of e* = 2. The 
dashed lines represent the analytical mean-field result from 
Eq. (3). The full widths at half maximum of both distribu- 
tions, denoted by AS and Ae respectively, are indicated by 
the thick gray arrows. 



potential Wshort('') = cT"^!^ jr and the mean-field inter- 
action by a Coulomb potential whose singularity at the 
origin is suppressed, W\av,^(r) = (1 — e~''/*)/r. For the 
screening parameter we choose s — ^RN^^^^ , i. e., much 
smaller than the initial nearest neighbour distance. With 
this choice VFiong closely matches the Coulombic case, 
while the initial interaction energy for Wshort is close to 
zero. Consequently, the inter-electronic repulsion energy 
resulting from Il^short cannot compensate the ionic back- 
ground potential Eq. (1) anymore. For realistically mod- 
eling with Il^short the effect of binary interactions on the 
PES, we drop the background potential but adjust the 
initial conditions such that the asymptotic single elec- 
tron energy of e* is preserved. 

The long-range case gives a final spectrum closely 
resembling the mean-field PES, albeit already slightly 
broadened at the edges. The short-range case, on the 
other hand, leads to a nearly symmetric spectrum sharply 
peaked at the single electron energy e* but with long tails 
reaching well beyond energies e = e* ± 1. Due to the 
rapid fall-off of the short-range potential most electrons 
do not interact with each other after photo-absorption. 
Only if the initial velocity vectors of two electrons put 
them onto a colliding trajectory an exchange of energy 
among these two electrons is achieved. Due to the high 
initial kinetic energy of the order e* only a small subset of 
initial conditions leads to electron pair trajectories with 
large energy exchange. Thereby, in a single binary colli- 
sion one electron can transfer all its kinetic energy to its 
collision partner and consequently the spectrum in Fig. 2 
covers the range e = . . . 4 for e* = 2. Since such violent 
binary collisions are very rare on the one hand, but lead 
to the largest energy exchange on the other hand, they 
can be viewed as an additional and largely independent 
random event, which augments the mean-field expansion 
dynamics. If this description is realistic, the convolu- 
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FIG. 3: Color online: (a) Final electron spectra for A*' — 10^ 
electrons propagated with Yukawa (blue, solid) and anti- 
Yukawa (red, dashed) potential, (b) Final spectrum from 
propagation with Coulomb potential (green, solid) and con- 
volution of Yukawa and anti- Yukawa spectra (black, dashed) 
according to Eq. (5). Dotted vertical lines indicate the excess 
energy e* = 2. 



tion Eq. (5) of the short-range and long-range spectrum 
should reproduce the full spectrum which is indeed the 
case as shown with the inset of Fig. 3. A residual inter- 
action effect manifests itself mostly in slight deviations 
at low energies: While the convolution of Eq. (5) ascribes 
each electron the same probability that its final energy 
e gets modified by a violent binary collisions, this is in 
reality more likely for slow electrons They come from the 
central, bulk-like region of the cluster and are more likely 
to suffer multiple collisions during their escape as com- 
pared to surface electron which are faster and suffer at 
most a single collision. 

The described construction of the PES from a mean- 
field and a binary-interaction component provides an in- 
tuitive physical picture. Yet, despite its approximate 
character it offers no computational advantage over the 
full result, since all trajectories for Wshort must be ob- 
tained numerically. Surprisingly, it is possible to take 
into account the correlation of mean-field and collision 
dynamics accurately by determining for each electron 
its pairwise isolated, binary Coulomb dynamics with all 
other electrons. This leads to the binary-intercation sum 
(BIS), a quasi-analytical and very accurate formulation 
for the PES which we introduce now. 

We consider for electron i the binary collision with elec- 
tron j, with initial positions and velocities (denoted with 
a prime) as in the Coulomb complex. Within the BIS 
approximation, the final energy of electron i is then 



energy Sij can be calculated analytically by means of the 
conserved Runge-Lenz vector in the center-of-mass frame 
[14]. Therefore we introduce the relative and center-of- 
mass coordinates r' = ^^'^ ^' — re- 
spectively. While the center-of-mass velocity is conserved 
(v = v'), we obtain an explicit expression for the final 
relative velocity v by means of the conserved angular 
momentum I = fi (r' x v') and the Runge-Lenz vector 
b = (v' X 1) + r'/r', with the reduced mass /i = 1/2. 
It reads v = — d (w (b x 1) — b) / (l -I- u^P) , whereby 
the absolute value v is known from energy conservation 

= + Finally, we use = {^r+^r/2f/2. 

The binary interaction does not explicitly include the 
influence of the background potential. However, it is easy 
to show that the definition (6) ensures conservation of the 
total energy which is, on the one hand, given by the 
l.h.s. of Eq. (6) E = J^i^i- ^^l^ other hand, it is 



e[.j) ~ — ^/''"ij from which follows that 



the r.h.s. of Eq. (6) is also equal to E = e'^ — V^ij- 
Note that the BIS is computationally extremely cheap 
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e[j and Sij are the initial 



with = £* by construction 

and final energies of electron i due the interaction or col 
lision with electron j, respectively, i.e., e - = ■i;^^/2-|-l/r^ 



with v[ = |v';| and r[j = |r^— r^| for two electrons with 
initial positions r^, and velocities v^, v^-. The final 



FIC. 4: Color online: Photo-electron spectra for propagated 
Coulomb complex (solid green/gray line) and BIS approxima- 
tion (dashed black) according to Eq. (6) for A*' = 10 (upper 
row), 10^ (middle) and 10^ (lower) on a linear (left column) 
and a logarithmic (right) scale, respectively. The gray ar- 
rows indicate the theoretical value of Ae = 3/4 according to 
Eq. (4) for the full width at half maximum of the mean-field 
distribution shown by thin solid lines. 
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FIG. 5: Color online: Photo-electron spectrum from a 2nm 
H2 cluster exposed to a pulse with 75 eV photon energy, cor- 
responding to E* = 60 eV (dotted vertical line). The pulse 
had an intensity of 2.5x lO^^W/cm^ and a duration of 500 as. 
The inset shows the same spectrum on logarithmic scale. 

since it does not require numerical propagation. Yet, 
quite a few analytical binary collision outcomes need to 
be summed since all N{N—1) pairwise interactions as 
well as multiple realizations of the isotropic velocity dis- 
tribution need to be taken into account. Hereby, it is 
crucial that the direction of an electron's velocity for all 
its A^— 1 binary interactions within one realization of BIS 
is kept fixed, to ensure the correct particle-particle cor- 
relations of the binary energy exchanges. 

The BIS spectra are shown in Fig. 4, overlayed with 
the fully propagated spectra. In all cases we see near 
perfect agreement. This agreement extends to specific 
features, such as a main peak arising from the residual 
discrete nature for the smaller systems with = 10 and 
N — 10^. Probably more surprisingly, also the long tails 
from violent collisions with very low probability are here 
reproduced in great detail (see the logarithmic spectra in 
the right column) . While these tails prevail for all system 
sizes shown, one sees a clear evolution towards mean-field 
dominated dynamics for larger systems with the central 
peak disappearing. In fact the shoulder on the left wing 
for the two smaller systems indicates roughly the maxi- 
mum due to mean-field dynamics. This is corroborated 
by the width of the distribution Ae, which is in good 
agreement with the value in Eq. (4) for all three cluster 



sizes if the shoulder of the spectrum is taken as the rele- 
vant maximum. 

For the sake of clarity we have restricted ourselves 
to possible single (photo-)ionization of each constituent 
of the cluster or large molecule and to sudden pho- 
toionization to introduce the phenomenon of massively 
parallel ionization. That this is nevertheless a realis- 
tic scenario is demonstrated in Fig. 5 with the PES for 
a 2nm hydrogen cluster induced by an XUV pulse of 
2.5x lO^^W/cm^ peak intensity, a duration of 500 as (fuU- 
width-at-half-maximum bandwidth of 4 eV) and a photon 
energy of 75 eV. This corresponds to an excess energy of 
E* — 60 eV (dotted vertical line in Fig. 5). Under these 
conditions, within reach by modern FEL sources, about 
10 % of the about 500 molecules in the cluster are ion- 
ized. The spectrum clearly shows the typical features for 
massively parallel ionization: a square-root shaped rise 
for energies E < E* and a high-energy tail for E > E*, 
see inset of Fig. 5. At very low energies one observes a 
structure due to electron-impact ionization which would 
increase for larger clusters. 

To summarize, we have introduced the phenomenon of 
massively parallel ionization which is a so far unexplored 
variant of multi-photon ionization and should routinely 
occur when illuminating larger targets with intense and 
short XUV to X-ray pulses. Characteristic for massively 
parallel ionization is a photo-electron spectrum which 
combines, almost independently, features from mean-field 
dynamics with those of violent binary electron collisions. 
Thereby, we could demonstrate that one can deduce from 
the width of the spectrum the ratio of charging (number 
of ionized electrons) to the size (radius) of the systems 
which is of great diagnostic value in experiments. More- 
over, we have devised a quasi-analytical yet highly ac- 
curate method to calculate the photo-electron spectrum 
from a sum of synthetic binary Coulomb collisions for 
small to large systems. Clearly, depending on photon 
energy and target, Auger-decay processes and multiple 
photo-ionization of one atom or small molecule within 
the cluster can occur and will modify the results pre- 
sented here. How the characteristic features of massively 
parallel ionization will be changed by such events will be 
investigated in future work. 
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